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We estimate the deformation of the cosmic microwave background radiation 

Ph| by the hot region ("cocoon") around a radio galaxy. A simple model by is 

O I adopted for cocoon evolution while the jet is on, and a model of evolution is 

'^ I constructed after the jet is off. It is found that at low redshift the phase after 

^ ' the jet is off is longer than the lifetime of the jets. The Compton y parameter 

generated by cocoons is calculated with a Press-Schechter number density 
rS ' evolution. The resultant value of y is of the same order as the COBE constraint. 

o3 \ The Sunyaev-Zel'dovich effect due to cocoons could therefore be a significant 

foreground source of small angular scale anisotropics in the cosmic microwave 

background radiation. 
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1. Introduction 

Ionized gas deforms the spectrum of the cosmic microwave background radiation 
(CMB) via inverse Compton scattering, the resulting spectral distortion being known as the 
Sunyaev-Zel'dovich effect (Sunyaev & Zel'dovich 1980; see for a recent review Birkenshaw, 
1998). The most popular source of the Sunyaev-Zel'dovich effect is hot X-ray emitting gas 
in a cluster of galaxies. The effects of the intracluster gas have been estimated theoretically 
by many authors (see e.g., Makino & Suto 1993;Aghanim et al. 1997), and there exist 
observations of the temperature decrement towards rich clusters of galaxies. Recently 
Aghanim et al.(1996) studied the effects of the partial ionization of the intergalactic medium 
(IGM) by quasar UV radiation. They showed that the kinematical effect due to the peculiar 
motions of the ionized bubbles is much larger than the thermal effect and can be detected 
by future generations of CMB observations. Natarajan & Sigurdsson (1999) exmanied the 
Sunyaev-Zel'dovich effect in the region formed by the quasar outflows. We examine here 
another source of the Sunyaev-Zel'dovich effect which is associated with radio galaxies. 

It has been suggested that a standard model for a typical strong double radio source 
comprises twin jets advancing into the intergalactic medium and surrounded by a hot region 
'cocoon' that consists of jet-supplied matter and shocked IGM (Scheuer 1974; Blandford 
& Rees 1974). Begelman & Cioffi (1989) and Nath (1995) constructed a model of the 
evolution of a cocoon in which the cocoon is overpressured against the IGM during its 
initial evolution. In their model, the expansion along the jet axis is determined by the 
balance of the thrust of the jet and the ram pressure, and the thermal pressure of the 
cocoon drives the expansion along the direction perpendicular to the jet axis. It has been 
demonstrated that some numerical simulations agree with this scenario (Loken et al. 1992; 
Cioffi & Blondin 1992). 

The energy density of the cocoon is higher than that of the intergalactic medium, and 
therefore it works as a local hot region that induces the thermal Sunyaev-Zel'dovich effect 
much as is observed for the hot intracluster gas in galaxy clusters. When we consider the 
lifetime of the jet(~ 10^ — 10® year), which is considered to be shorter than the Hubble 
time at the redshift where radio galaxies reside, cocoons may not play a significant role in 
the deformation of the CMB spectrum. However, even after the jet energy supply stops, 
cocoons can remain hot; the cocoon expands due to its own thermal pressure, and the 
thermal energy of the cocoon is consumed by expansion and radiative cooling by the cocoon 
gas. Thus the cocoon remains a possible source of the Sunyaev-Zel'dovich effect for an 
effective lifetime that can be much longer than the jet lifetime. 

We estimate the Sunyaev-Zel'dovich effect due to the cocoons by taking into account 
their evolution both while the jet is on and after the jet is offi We make use of the simple 
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model of Nath (1995) for the evolution while the jet is on, and first estimate a rough 
order of y for a single cocoon. We construct a model for the evolution after the jet is off, 
and evaluate the average Compton y parameter for a uniform distribution of cocoons. 
Comparing our result with the COBE constraint, we discuss the possibility of detection of 
the effect of cocoons in the CMB observations. 

Our paper is organized as follows. In §2 we summarize our basic calculation and list 
the parameters we need to estimate in order to apply evolutionary models for the cocoons. 
In §3 and §4 we describe models of cocoon evolution during the period while, and after, the 
jet is on. In §5 we calculate the Sunyaev-Zel'dovich effect and give our results. Discussion 
and conclusions are given in §6. 
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2. Basics and Necessities 



First we represent the formal expression of tlie Compton y parameter and clarify what 
parameters we need to estimate. For a single cocoon, 



kT 



f kl 
l/sc = / :TneaTdR, 

^-tV''\ (1] 



p. 



rripC^ 



where o"t, Ug, T, Pg, dR are the Thomson scattering cross section, electron number density, 
electron temperature, pressure, and distance element along the line of sight, respectively. In 
the second equation, we set the length scale along the line of sight R ~ V^^'^ as averaged for 
directions over a randomly oriented distribution of cocoons. The total y in a line of sight 
with a solid angle element dVL is the integral of the product of |/sc, the formation rate of 
cocoons, the volume element a^r'^drdQ, and the crossing probability, as 



yUM, t - tcoi)^(M, tcoi)dtcoia\^drdn—dM, 



/ D \ 2 

ysc{M, z)n{M, z)a^r^dr ( — ) dM, 



PV 1 

^—aTniM,z)a^r'^dr—^dM. (2) 

rrieC^ Ri 



Here, a is the cosmological scale factor, r is the comoving radial coordinate, Ra is the 
angular diameter distance, 6 = R/ Ra-, n is the number density of radio galaxies with mass 
M at tcoih d"^ /dVt is the probability for the line of sight to traverse the cocoon, respectively. 
The angular diameter distance Ra is given by 

Ra = a{z)r = jj^^[l - (1 + z)-'/']. (3) 

[i + Z) Uq 

Strictly speaking, other than their mass scales, the epoch of the collapse of objects depends 
on various factors such as the amplitude or density distribution within each object. 
However, in the second and third equations, we assume for simplicity that the objects 
which have the same mass collapse at almost the same epoch, and replace dn/dt x dtcoi 
with n{M, z). Thus in order to evaluate y, we only need the value of PeV for each cocoon. 
In §3 and 4 we examine the evolution of PgV within a cocoon with the jet on and off. 
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3. Model 1: while the Jet is On 

In this section we consider an evolutionary model of a strong double radio source with 
jets (Scheuer, 1974; Blandford & Rees 1974; Begelman & Cioffi 1989; Nath 1995) and 
make a rough estimate of y for a single cocoon. We follow the model proposed by Nath 
(1995), in which the jets are enveloped in a hot cocoon consisting of shock-compressed 
jet matter and intergalactic medium (IGM). In this model, the jets are confined by the 
pressure of the cocoon, i.e., we do not take confinement by magnetic pressure into account. 
Highly collimated, 'light' (the density of the jet matter is smaller than the ambient IGM) 
jets interact with the IGM and form bow shocks. At the shock, some of the electrons are 
accelerated to high energy and radiate synchrotron emission. The thermalization efficiency 
is one of the unknown parameters in this model; however, because in many lobes of radio 
galaxies, magnetic fields have been observed, we expect that the effective mean free path 
of the particles should be short and the thermalization efficiency should be high because 
of the interaction with the magnetic field. Therefore in this paper, the relativistic jet 
matter is assumed to be completely thermalized at the interaction interface with the IGM. 
The thermalized jet matter forms a high temperature region around the jet together with 
the shock-compressed intergalactic medium. The shock-compressed intergalactic medium 
surrounds the thermalized jet matter separated by a contact discontinuity (see Figure |l|). 
At least initially, the thermal pressure of this hot region is higher than the IGM pressure. 
When the hot region pressure is higher than that of the IGM, the hot region also expands 
in the direction perpendicular to the jet axis and is accompanied by shocks. We make use 
of the term 'cocoon' for the hot region surrounded by the shock surface. As the cocoon 
expands the pressure decreases. Nath (1995) estimated the equiUbrium timescale, i.e., the 
time required until the cocoon pressure decreases sufficiently with time to be equal to the 
ambient pressure: 



t.„ ~ 1.34 X 10^° r^^ ^V'lfhl"( ^ ^ .^l^i^h 




•^^ — -^" Vi06K; "IGM -50 l^io45ergs/sec; '^^' ^(SOkpc)^; \l/'i 

(4) 
where Lj is the kinetic power of the jet, Vh is the shock advance speed into the IGM, 
j3j = Vj/c is the jet speed. Ah is the effective cross section of the bow shock, e„ is the volume 
factor (which depends on the shape of the cocoon), ^^igm is the density of the IGM in unit 
of the critical density, Hq = 50/i5okm/sec/Mpc, respectively, and the suffix a denotes IGM 
value. It is found with observed values of Cygnus A that teq is longer than the jet lifetime 
as long as Ta ^ lO'^K: i.e., cocoons remain overpressured while the jet is on. 

Hereafter we consider two collimated steady jets advancing into the ambient IGM. 
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Fig. 1. — A schematic diagram of our model of the cocoon surrounding the jet while the jet is on. See text 
for details. 

The thermalized jet matter and the shock-compressed IGM matter form a cocoon around 
the jets and the cocoon expands with shocks advancing in directions both parallel and 
perpendicular to the jet axis. Below we will derive some relations between the various 
variables in order to obtain PcV during this stage of evolution. 

At the shock, the jet thrust is balanced by the ram pressure of the IGM at the shock 
that works on the effective area A^, 



AhPavl 



(5) 



The kinetic power Lj is written as follows with jet speed Vj. From observations and theory 
of the superluminal motions of the quasars, Vj is considered to be approximately equal to 
the speed of light. Therefore, 



L, 



^Pjv^^j^ 



(6) 
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Vj = c. (7) 

Here Aj is the cross section of the jet. At the interaction of the jet and IGM, a bow shock 
is formed. Following Nath (1995), we approximate the bow shock as two oblique shocks. 
Then the jump condition at the oblique shocks is 

Pc 5{vh/cafsm^(j)-l 

y. = 1 ' ^'^ 



where is the angle between two oblique shocks, and Ca is the sound velocity of the ambient 
medium. 

As we mentioned above, because the cocoon is overpressured against the IGM, the 
shock also advances into the IGM in the direction perpendicular to the jet. The jump 
condition there is 



where Vc = dr^/dt is the shock advance speed. 



l-r'r^PaV^, (9) 



Finally, we assume that the internal energy within the cocoon is supplied by the jet 
matter, and that jet kinetic energy is transformed into the thermal energy of the cocoon gas 
and the kinetic energy to expand the cocoon. We also assumed the fraction of input energy 
that is converted to other forms of energy (magnetic fields or high energy cosmic rays: see 
discussion in §6) can be neglected. From strong shock conditions, Vc is the same order as 
the random velocity of the cocoon gas, and therefore the fraction of the thermal energy to 
the total input energy t/iGM/-E'mput is 0{1), 

-^PcV ^ L^t, (10) 

7-1 

V = ey.Vht, (11) 

where 7 is the adiabatic index and t is the time elapsed since the jet is on. The final 
equation comes from the assumption that Vh is constant Ih ~ Vht. Numerical simulations 
indicate that Vh decreases and Ah increases with time (Loken et al. 1992; Ciofii & Blondin, 
1992). However, the time-dependence is relatively small (//, oc t" with a ~ 0.8 in their 
simulations), and therefore we assume Vh is constant for the estimate of ysc- Note that for 
the estimate of total y, as described in §2, we only need PcV, and therefore do not need the 
detailed evolutionary scenario (see eq.[|lO|). 



In order to estimate the Compton y parameter value, we need the number density 
n{M, z) and PcV of objects with mass M and redshift z. As we show in §5, we use the 
Press-Schechter formalism for n{M,z). Therefore we need to write the product of the 



pressure of the cocoon P^. and V as function of t, mass of the host galaxy M, and the 
parameters of the ambient medium. We make additional assumptions about the jet kinetic 
luminosity. We assume that Lj is equal to the Eddington luminosity of the black hole 
-^edd = 4:7rcGMBK/Ke ^ 1.25 X 10^^ {Mbr/Mq) ergs/sec. For the mass of the black hole 
some recent observations showed that there is a rough proportionality between the mass 
of the black hole and the host galaxy, Mbh = 0.002M (Magorrian et al. 1998; Ford et al. 
1997; van der Marel 1997). Using these assumptions we can write Lj in terms of the mass 
of the host galaxy. Integrating equation (||), we obtain the relation 



r? = 2 f I^ A. 



1/2 



t . 



\ (^v PaVhJ 

Using the above equations we obtain the following expressions for P^ and V, 
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(13) 
(14) 



According to the usual discussions of synchrotron spectral aging, the minimum age of 
the jet is typically about 10^^'^ years (Liu, Pooley, & Riley 1992). So if we take cosmological 
values for the ambient medium with Qhh"^ = 0.01, 



Pc ~ 2.13 X 10 
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Mpc, 
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and y, 
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' '-^ Y'\l + zfl^ l^^\''\ (17) 



Thus we expect that the deformation of CMB caused by cocoons is possibly relevant to 
observations. 
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4. Model 2: after the Jet is Off 

After the relatively short lifetime ~ 10^ years, the jet turns off. However, after the jet 
is off, cocoons may retain a high temperature and remain a source of the Sunyaev-Zel'dovich 
effect for a long time. A cocoon cools radiatively and expands due to its own thermal 
overpressure. We need to estimate the effective lifetime of the cocoon with the jets off, as a 
'relic' of the radio galaxy. 

A cocoon can be considered as an energy input into the intergalactic medium. Thus 
the evolution of a cocoon after the jet is off resembles the evolution of a supernova remnant 
in the interstellar medium. We may therefore adopt the analogy with the evolution of a 
supernova remnant. It is a well-known result that there are three stages in the evolution of 
a supernova remnant in the interstellar medium (Spitzer 1978). First, the remnant expands 
freely at nearly constant speed into the surrounding medium. In this phase, the deceleration 
due to the swept-up interstellar matter in front of the shock does not intervene in the 
expansion. As the shock plows into the interstellar matter, it begins to decelerate. When 
the mass of the swept-up matter (47r/3)i?^pisM becomes nearly equal to the mass of the 
ejecta Mejccta, one needs to take into account the deceleration of the shock front. Thereafter 
the evolution goes into the next stage, the Sedov (adiabatic) phase. In this Sedov phase, the 
temperature behind the shock front is so high that deceleration due to cooling behind the 
shock front can be ignored compared with the expansion. There is a well-known similarity 
solution for the evolution of the shock front. 



Rs = UEt^lpiSMf'\ (18) 

Vs = 0.4eo(^/pisMt')'/', (19) 

where E is the initial energy input and .^o = 1-17 with 7=5/3 (Shu 1992). As the shock 
decelerates and the temperature behind the front decreases, cooling begins to affect the 
expansion and a dense shell is formed behind the shock front. Almost all the mass of the 
ejecta and plowed-up interstellar medium are in the shell, and the thermal energy of the 
remnant is radiated away by the shell. The condition to be satisfied at the transition from 
the Sedov phase to the radiative phase is Tex = Rs/vg ~ Tcooilshock- When this condition is 
satisfied, i.e., the cooling time subsequently becomes shorter than expansion time of the 
shell, the thermal energy of the remnant is radiated very rapidly (Chevalier 1974). 

We will examine the evolution of the cocoon after the jet turns off by analogy with 
the evolution of the supernova remnant described above. Cocoons should have three 
evolutionary stages as do supernova remnants. Let the cocoon evolve into the Sedov phase 
at t = tg and into the radiative phase at t = t^ with t=0 at the time when the jet stops. 
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First, the cocoon expands with nearly constant speed. The shock surrounding the cocoon 
engulfs the intergalactic medium. Accordingly the kinetic energy is turned into the thermal 
energy, and the fraction of the thermal energy of the shock- compressed IGM to the total 
energy input by the jet, f/iGM(^)/(-Z^j^iife), increases with time. Let us consider similarly as 
the case of supernova remnant further. The deceleration by the piled-up IGM cannot be 
ignored when the fraction of kinetic energy, \AEk/Ek\, which is transformed into thermal 
energy, approaches unity. At this epoch we can consider that the Sedov phase sets in. As we 
wrote in §3, while the jet is on, t/icM is the same order as Ek, and therefore one can regard 
that the Sedov phase begins when \/S.Ek{ts) / Ek{Q)\ ~ At/iGM(^s)/f^iGM(0) = Oil). In other 
words, ts is about the same as the time required to double the thermal energy of the cocoon 
at the time when the jet stops (t = 0). Remembering that the cocoon expands almost 
constantly before the Sedov phase because of lack of swept material which decelerates the 
expansion, i.e., A?7iGM(i)/At ~ Lj ~ f/iGM(0)/tiife, we can consider that the cocoon evolves 
into the Sedov phase at t = t^ ^ tufe, since Af/iGM(^s) ~ f^iGM(O). 

The radiative phase begins when the cooling time behind the shock front Xcooilshock 
catches up with the expansion time of the Sedov phase Rg/vg. Because the thermal energy 
of the cocoon is radiated away very rapidly once cooling sets in, we expect that the cocoon 
remains hot enough to be the source of the Sunyaev-Zel'dovich effect only during the first 
two phases. Therefore in effect the lifetime of the cocoon can be taken as ~ t,.- 

In order to know U, we need to evaluate Xcooil shock and hence we need to know the 
cooling processes of the cocoon. First we examine the content of the cocoon after the 
jet is off. The total thermal energy of the cocoon is the sum of that of the jet-supplied 
matter and of the shock-compressed intergalactic matter. Within the model described in 
the previous sections, the jet-supplied matter and the compressed intergalactic matter are 
in pressure equilibrium, and we therefore cannot distinguish the former from the latter. 
Thus we do not know the ratio of the thermal energy of the jet-supplied matter to the 
total energy. Hereafter we assume that the thermal energy of the shock- compressed matter 
already dominates that of the jet-supplied matter at the time when the jet stops. 

We need to include various cooling processes that affect the temperature behind the 
shock front, 

Te = ^. (20) 

From the discussion of synchrotron spectral aging, the shock speed Vc is estimated to be 
Vc < Vh ^ 0.1c (Liu et al. 1992). This result means that T^ is lower than ~ 10^° K from 
equation (^). We can thus ignore pair creation (e~ — e"*") within the cocoon. We regard 
the cocoon matter as a hydrogen plasma {p + e^). The dominant cooling processes within 
this plasma are thermal bremsstrahlung and inverse Compton cooling, and hereafter we 
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Fig. 2. — Effective lifetime of the cocoon after the jet is off tr against redshift. At higher redshift Compton 
cooling is more effective and tr is shorter. 

take account of these processes for the estimate of TcooiI shock- 

In figure 0, we plot t^ versus redshift. At higher redshift, inverse Compton coohng is 
more efficient and the coohng time is shorter. As noted above, ts is about the same order as 
the jet hfetime tufe ~ 10^ — 10*^ years. This means that when tr < tufe, tr < ts and no Sedov 
phase would exist. In this case we set tr = 10^ years (the final result is hardly affected by 
this value). On the other hand, when tr > tufe, tr > ts and the Sedov phase exists. We make 
a fitting formula for tr as a function of redshift: 



< 2; < 5 

5 < z < 10 

10 < 2 < 20 



logio(i^r/years) = 10.7 — 0.37 z , 
\ogiQ{tr/ years) = 9.93 — 0.21z , 
logio(^r/years) = 8.96 — O.llz . 



The thermal energy of the cocoon is expected to decrease rapidly at t > t^ so we assume 
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P,V^exp(^-^y (21) 



after the jet stops. 



5. Averaged y Estimate 

In the previous sections, we have obtained all the building blocks for a single 
cocoon in order to calculate the total y. One remaining parameter is the cocoon 
number density n{M,z). We employ the Press-Schechter formalism, assuming 
n{M,z) = fr{M,z) X nps{M,z), where fr{M,z) is the fraction of radio galaxies 
among the objects which have mass M at z. In the next subsection we will model fr{M). 



5.1. Estimation of fr 

As we described in §4, cocoons can remain hot for a long time. Therefore in order 
to calculate y, we should count the "relics" of radio galaxies. However, we do not have 
any theoretical model of fr{M): instead we make a rough estimate of fr{M) from the 
present-day luminosity functions of radio galaxies and field galaxies.0 

First, for radio galaxies, we make use of the local luminosity function of radio sources 
of Condon (1989). His sample contains both the active galaxies with massive black hole 
( "monsters" ) and starburst galaxies. But these two types of galaxies cover quite a different 
range of radio power: most of the extragalactic radio sources with P(1.4GHz) < lO^^W 
Hz~^ are starburst galaxies and most of the sources with P(1.4GHz) > lO^^W Hz~^ are 
"monsters" (see his figure 9; see also Sopp & Alexander 1991). 

Next, in order to compare the numbers of radio galaxies with ordinary galaxies, we 
examine the (optical luminosity)- (radio power) relation for radio galaxies. We employ the 
data of Zirbel & Baum (1995) and plot the P(2.7GHz)-absolute V magnitude relation.^ 
From figure ^, we can see that there is an upper limit of absolute magnitude My ~ —22 
mag though there is a large scatter. This means that there is a small number of radio 
galaxies fainter than —22 mag. Finally, for the luminosity function of field galaxies we 



^In this subsection we use h = 0.5 and qo = 0.0 in order to use the resuhs in the Hterature. 

■^We make use of the power law index a in the Table of Zirbel & Baum (1995) to calculate P(2.7GHz) 
from P(408MHz), and when a is not displayed in the Table, we set a = 0.75. 
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Fig. 3. — Radio power at 2.7 GHz and absolute V magnitude relation for radio galaxies. Data is taken from 
the table of Zirbel & Baum (1995). Open squares are FRl galaxies, filled squares are FRII sources, crosses 
are intermediate ones, triangles are those with unidentified morphology, and open stars are high z(> 0.5) 
galaxies. Zirbel & Baum (1995) adopted h = 0.5 and go = 0.0. 

adopt the Schechter function. Note that usually the Schechter function is a function of Mb, 
not of My- Therefore we have to obtain the My-Schechter function. We include 4 Hubble 
types of galaxies, E/SO, Sab, Scd,and Sdm, and adopt the luminosity-averaged local fraction 
of each types from Yosii & Takahara (1988). We adopt the value ifj* = 2.3 x 10^^ Mpc~^, 
M^ = —21.1, and a = —1.11 from Yosii & Takahara (1988), and the averaged B — V colors 
for each morphological type of galaxies are taken from Babul & Ferguson (1996). The 
resulting My-Schechter function is displayed in figure |[ Therefore we obtain a model of 
the fraction of radio galaxies in terms of 



luminosity, 



friMv < -22mag) 
fr{Mv > -22mag) 



/ 0(-Pradio)'^-Pradio 



0, 
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''Schechter 



{My)dMy 



0.01, 



(22) 
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Fig. 4. — Schechter function in terms of My and Mb- 

where is a luminosity function of radio galaxies and 0schechter is the My-Schechter function. 

Because we use the mass function of galaxies to count the cocoons, we need /,. as a 
function of mass. From our observation of figure ^ that radio galaxies have a lower bound 
to their optical luminosity, we further assume that radio galaxies have a lower mass limit 
corresponding to My ~ —22 mag. The critical magnitude My = —22 mag is about the 
same as Ly {B — V^ 1-0), and the lower mass limit to the total (baryons+dark matter) 
mass of a radio galaxy can be taken to be ~ IO^^Mq. It is natural to consider that fr 
depends only on masses greater than IO^^Mq but moreover we assume that for simplicity 
fr is independent of M in this mass range. Then we can model fr{M) as 



friM) = /; = friMy < -22mag) = 0.01 M > IO^^Mq, 

fr{M) = M < lO^^M©. 



(23) 



From the popular idea that the host galaxy of a radio galaxy is a giant elliptical, we 
set the upper limit on the mass to be IO^^Mq. 
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5.2. Integration and Results 



The physical number density of cocoons is proportional to 



nps{M,z)dM 



^p{z) S, 



n M Diz) 



1 da 



a- 



dM 



exp 



2a^D{z)'^ 



dM, 



(24) 



where p{z) is the proper density at z, 5c = 1.69, D{z) is the growth of the perturbation 
{D[z) = (1 + z)~^ for the Einstein-de Sitter universe), and a^ is the variance with a 
Gaussian filter. 



a\M) 



1 r°° 

iTT^ JO 



(25) 



Here tm = {M/ pq)^^^ /^/27r is the smoothing radius of the Gaussian filter with present 
critical density po, and P{k) is the power spectrum of density perturbations. In calculating 
P{k) we use the fitting formula of Bardeen et al. (1986), 



P{k) = AT{kfk, 

T{k) = Ml-tA^[i + 3.89g+(l6.1g)2 + (5.46g)3 + (6.71g)T'/', 



2.34g 
k 



(26) 

(27) 



Q 



nh'^Mpc 



-1 ■ 



(28) 



In this paper we set Q = 1.0 and h = 0.8. For the normalization factor A we use the 
CO-Bi?-normalized value (Sugiyama 1995) and approximate 



1 da 1 [{n + Sy 



a^dM M\ 6 



(29) 



with P{k) oc A;" within the range of mass of a radio galaxy. With IO^'^^Mq < M < IO^^Mq, 
—2.5 ^ n ^ —0.5 so we set n = — 2 throughout. 



Thus the total averaged y is given by 

kT 



y 



single cocoon 171^0 



nardR frnps{M, z 



\ ' ^ a^r^drdMdn, 
dil 



PrR 



nipC^ 



arfrnpsiM, z)e\M, z)a:h-'drdM. 



(30) 
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Therefore, 



-") 



f P W 1 

y = fr^-^crTnpsiM,z)-^a^r^drdM, 

If we use equations (|T^ and (0) , 

PV={^^^^'~^^^' :^<%e .32X 

\ Lj(7 - l)tHfe X exp(-(t - tiifc)/t,.(2;)), :t > %« 



applying equation (|21|) . 



We integrate equation (|3T|) with lO^^M© < M < 10^'' M© and < ;z < 20. For the 



hfetime of the jet we take tufc = 10''' years regardless of mass. We obtain 



,_5 ( /bh \ / fr 
0.002/ I 10-2 



I/tot ^ 5.72 X 10-^ -^ kfe • (33) 



As we discussed in §4, cocoons serve as the source of the Sunyaev-Zel'dovich effect after the 
jet is off, and the number of the cocoons can be greater than the number of the luminous 
radio galaxies. 

Our estimate of y should be compared with the COiJi? FIRAS constraint y < 1.5 x 10~^ 
(Fixen et al. 1996). Discussions about our result and this constraint are presented in the 
final section. 
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6. Discussion 

We have investigated tlie Sunyaev-Zel'dovicli effect produced by cocoons of radio 
galaxies. We liave employed a simple model of Nath (1995) while the jet is on. We have 
also examined the evolution of a cocoon after the jet turns off, adopting the analogy to 
the evolution of a supernova remnant in the interstellar medium. The cocoon remains 
at high temperature after the jet turns off for over ^ 10^ years to low redshift z <^ 10, 
remaining overpressured against the pressure of the IGM. Because of the short lifetime of 
the jet (~ 10^ — 10^ years), the main contribution to y comes from the cocoons with the 
jet off. Comparing our result with the COBE constraint we find that hot cocoons could be 
significant foreground sources of the CMB. We discuss the assumptions and the implications 
of the models below. 

We have made a model for fr{M) which is a step function with a cut-off on the upper 
end of the mass range. We assumed the lower limit on the mass of a radio galaxy to be 
IO^^Mq from the optical magnitude distribution for radio galaxies, but this value has some 
ambiguity. We performed calculations changing the lower mass limit Miow, and found that 
y increases very rapidly with decreasing Miow to exceed the COBE constraint by far. This 
might support the view that there exists a low mass limit for jet activity. 

In §4 we have presumed that the cocoon remains overpressured against the intergalactic 
medium. However, we must investigate the possibility that the pressure of the cocoon 
decreases to be in pressure equilibrium with the surrounding medium within the effective 
lifetime of the cocoon ~ tr- Once in pressure balance with the IGM, the internal energy 
of the cocoon is consumed by PdV work. We calculate the timescale to be in pressure 
equilibrium with IGM tp as a function of redshift. There has not been any direct 
observational means for obtaining the pressure of the IGM. As for the temperature Tq, one 
can expect T^ ^ 10^ — 10^ K from the previous studies of the reionization of the Universe 
(see e.g., Fukugita & Kawasaki 1994;Tegmark, Silk, Sz Blanchard 1994). However, radio 
galaxies or radio loud quasars are likely to be found in high density regions (Jones et al. 
1997). Pressure in such an environment is higher than in an ordinary region. Thus we take 
the temperature of the intergalactic medium Tiqm to be 10^ — 10^ K, on the order of the 
present value of the temperature of the intracluster gas. For the number density of IGM 
Ua, as an extreme case, we equate n^ with the critical density ricr = ?>Hl / {^TiGmp){l + 2;)^, 
which is much larger than the average baryon density by a factor of Vtbh"^. On the other 
hand, we make use of the similarity solution for the Sedov phase to calculate the pressure 
behind the shock front. In figure ^ we plot tp against redshift z. From this figure we 
see that the lifetime of a cocoon in a rich cluster of galaxy is shorter than tr- For the 
most extreme case, with the electron number density rig ~ 10~^cm~^ and the temperature 
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Fig. 5. — Timescale to reach the state of pressure equihbrium with the external IGM pressure tp against 
redshift. The temperature of the IGM is constant. Sohd hne corresponds to the case with Tiqm = 10^ K 
and dashed hne corresponds to the case with Tjgm = 10^ K. 



T ~ lO^K, tp is about 1 x lO^years. We find that when U > lO^years for z < 10, 
y ;^ 2.6 X 10-^(/bh/0.002)(/^/10-2), which is still observable. 

For radio luminous galaxies, the radio emission itself may affect the Sunyaev-Zel'dovich 
decrement. In order to estimate this effect we should take into account the integrated 
spectrum in a line of sight from cocoons at different redshift. From equations (17) and 
(33), however, we see that there is marginally one cocoon in one line of sight, which means 
that the contamination by radio emission from other radio galaxies is not a severe problem. 
Thus we only have to consider the radio emission accompanying the cocoon itself. Detailed 
spectrum evolution requires a specific model of radio lobes, though, we just present simple 
discussions in this paper. 

The Sunyaev-Zel'dovich decrement of fiux is (Zel'dovich & Sunyaev, 1969) 



AL 



xe 
fe- - 1] 



X 



■V: 



(34) 
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with the nondimensional frequency x ~ hv/kT . In the Rayleigh- Jeans region(x -C 1^ 

A/. 



/, 



-2y, (35) 



V 



which is independent of x, and the CMB flux containing Sunyaev-Zel'dovich signal has also 
black body spectrum. Kt v '^ lOGHz, the spectrum of radio galaxy is I^, ex t/~" (typically 
a > for bright sources; see Sokasian et al. 1998), which is quite different from the CMB 
spectrum (oc z/^). As Sokasian et al. (1998) pointed out, bight radio galaxies can be major 
sources of the anisotropy of CMB especially on subdegree scales at lOGHz^ z/ ^200GHz, 
and the Sunyaev-Zel'dovich singal can be completely overcome by the radio emission 
in this waveband and angular scale while lobes are bright. In the Wien region, radio 
emission from radio galaxies becomes increasingly fainter with v, and the influence on the 
deformation of CMB is far weaker than in the Rayleigh- Jeans region. Thus, one coule 
consider that the observational prediction of the integrated y of the Sunyaev-Zel'dovich 
effect is significantly contaminated by emission from the lobes in low frequencies. In such a 
case, the multiwavelength observation covering a wade range including both Rayleigh- Jeans 
and Wien region, which PLANCK satellite will carry out^, is required in order to detect the 
Sunyaev-Zel'dovich effect. 

However, the effect of radio emission from the lobes depends on the ratio of the 
timescale of energy loss of the high energy electrons by synchrotron emission and the 
cocoon effective lifetime(~ tufc + t^). Synchrotron emission from lobes may not fade away 
immediately after the jet stops, while cocoons can survive for a relatively long period. If the 
lifetime of synchrotron emission is significantly shorter than the effective cocoon lifetime, we 
expect little contamination from radio emission to measure the Sunyaev-Zel'dovich signal 
in the Rayleigh- Jeans region. The loss timescale can be estimated if we know the magnetic 
field strength and the density of high energy electrons at the lobes. Note that McKinnon, 
Owen, & Eilek (1991) reported the detection of Sunyaev-Zel'dovich decrement from radio 
jet lobes with NRAO 12m telescope. 

The lobe decay time is also closely related to the fraction of high energy electrons e* 
(see below) . However, the estimate of the lobe decay time requires a detailed model of radio 
emitting lobes, and this is beyond the scope of this paper. In addition to the evolution of 
the spectrum, the estimate of lobe emission decay time should be regarded as a future work. 

Radio galaxies have one or two lobes that emit synchrotron radiation from high energy 
electrons, which are accelerated via Fermi acceleration at the shock. We assumed in the 



see http://astro.estec.esa.nl/SA-general/Projects/Planck/ 
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cocoon model that almost all of the jet matter is completely thermalized at the shock, but 
more reahstically speaking, only a fraction ej < 1 of the jet-supplied energy can become 
the thermal pressure of the cocoon. Therefore Lj in equation (|52|) is replaced by ejLj and 
y is ety. As expressed in §3, when the jet is stopped by the shock, a large fraction of jet 
matter would be thermalized (et ~ 1) with the magnetic field. The thermalization efficiency 
et can be estimated if we know the species of the jet particles (whether they are neutral 
particles or electron-positron plasma as expected by the central engine of a massive black 
hole) and the acceleration mechanism at the shock. X-ray emission from the lobes, which 
is considered to be inverse Compton scattered synchrotron photons, will provide important 
information about Cf. 

Jones et al. (1997) found the decrement of CMB temperature towards the z = 3.8 
radio-quiet quasar pair PC1643-I-4631 A and B. There is no detected optical or X-ray 
source in this direction, and they consider that a low surface brightness cluster of 
galaxies may be present at relatively high {1 < z < 3.8) redshift. Another explanation 
is suggested by Natarajan & Sigurdsson (1999), who consider that the decrement is the 
kinematical Sunyaev-Zel'dovich effect produced by quasar winds. Neither explanations 
seems particularly plausible. We suggest that the so-called "black cloud" phenomenon (see 
also Richards et al. 1997) may be related to the relics of radio galaxies, with the lobes 
already vanished and with the cocoon still hot enough to give a sizable Sunyaev-Zel'dovich 
decrement. 

We thank H. Susa and R. Nishi for useful discussions in the early stage of this study, 
and anonymous referee for helpful comments. M.Y. is grateful for many valuable comments 
by K. Ohta about the luminosity functions. The research of J.S. has been supported in 
part by a grant from NASA. This work is financially supported in part by the Grant-in-Aid 
for Encouragement of Young Scientists by the Ministry of Education, Culture, Science and 
Sports, No.09440106. 
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